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ABSTRACT 

Crop identification using multispectral satellite imagery and 
multivariate pattern recognition is a relatively new technique enabling 
rapid evaluation of large areas. The accuracy of this process is 
examined in this study in a semi-arid climate. Multispectral reflec- 
tance data was collected by the Earth Resources Technology Satellite 
(ERTS-1) over the semi-arid regions of western Kansas and Te'-as during 
the 1973 growing season. 

Multivariate pattern recognition was used to identify wheat 
accurately in Greeley County, Kansas. A classification accuracy of 97% 
was found for wheat and the wheat estimate in hectares was within 5% of 
the USDA’s Statistical Reporting Service estimate for 1973. 

The multispectral response of cotton and sorghum in Texas was noc 
unique enough to distinguish between them nor to separate them fro^-A 
other cultivated crops, either singly or multitemporally . The teat 
site of Lubbock County, Texas was deemed too heterogeneous in agricul- 
ture practices for correct identification of cotton and sorghum using 
ERTS-1 imagery. 

ERTS-1 Imagery may be a useful tool In improving crop surveys. 
Current data acquisition systems and analysis techniques worked quite 
well for a homogeneous agricultural area like Greeley County, Kansas. 
Areas which are quite heterogeneous in agricultural practices, crops 
and soils are problem areas for which acceptable crop identification 
may not be obtainable using ERTS-1 imagery. Area estimation of crops 
In heterogeneous regions does not seem feasible using present satellite 
imagery . 


CHAPTER I 


INTRODnCTION 


Crop surveys have been considered veary lo^ortant in Che past few 
decades. Their importance will probably remain quite relevant as our 
total population Increases. Current, accurate crop surveys could help 
stabilize supply-demand relationships for farm products. Distribution 
of produce from agricultural crops would be more timely If the concen- 
tration of crop production was known. The producer , the processor and 
the couBumer of crop products would all benefit from accurate crop 
surveys over the long run (Eisgruber, 1972) . The chaos created by food 
shortages would be lessened. 

Various federal agencies, state agencies and private organizations 
have strived to obtain crop surveys which would help in forcastlng the 
crop production of a given area. Most surveys to date have been based 
on information volunteered by farmers and ranchers and on a very small, 
random sample of the crop land which Is observed in the field and 
mea. lured by field personnel. The accuracy of such crop surveys In the 
pa&t has shown that crop surveying methods could be Improved (Gumielson, 
et al., 1972). However, to increase reliability using the crop survey 
techniques of the past, the number of field observations should be 
increased; a very costly undertaking. In this thesis a method of 
sampling much larger areas of crop land at reduced costs will be 
examined; specifically identification of agricultural crops with 
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multlspectral satellite iioagery using multivariate pattern recognition 
techniques. 

Remote multlspectral sensing (Lars, Vol. 3.) may be defined as the 
sending, from: a remote location, of electromagnetic radiation - either 
reflected or emitted - in many discrete, usually relatively narrow 
spectral bands between wavelengths of 0.3 pm and 15 pm and also in radar 
bands from about 0.86 to 3.0 centimeters. 

The new view of the Earth from the Earth Resources Technology 
Satellite (ERTS-1) at 940 kilometers altitude may provide timely and 
accurate information which could lead to more useful crop surveys. 

Since its launch in July 1972, BRTS~1 (NASA hoe renamed BRTS-l; now 
LANDSAT-1) has been scanning the entire Earth every 18 days. The 
satellite views a swath 190 kilometers wide as it orbits the earth. 

The spacecraft carries a multlspectral scanner to obtain image data 
in various spectral ranges (green, red and infrared). More than 5000 
images covering about 180 million square kilometers are collected each 
week. Techniques to analyze ERTS-1 imagery by computer have been 
developed at the Laboratory for Applications of Remote Sensing (LARS), 
Purdue University and these techniques will be discussed further into 
the text. It was felt that the wide-area, sequential coverage of the 
ERTS-1 imagery, combined with the capabilities of computer processing, 
offered a new opportunity to identify crops and Improve crop surveys 
over large areas. 

Two areas were selected in the southern Great Plains as trial 
sites; the sites being Greeley Cotmty, Kansas and Lubbock County, Texas. 
Both sites are considered semi-arld. The major crop of Greeley County 
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Is winter wheat and the loaln crops of Lubbock County are cotton and 
grain sorghum. If these crops could be correctly Identified and mapped 
using ERTS-1 Imagery, which samples all fields in a given area, then 
this information should demonstrate the usefulness of ERTS-1 imagery. 

A more complete sampling of an area should result in more accurate 
crop estimates than are currently available from conventional surveys. 

In summary, the overall objective was to determine if crops could 
be correctly Identified using ERTS-i imagery . A complete, reasonably 
correct Identification of crops for an area should provide additional 
useful Information which should make possible more accurate estimation 
of crop production. 


CHAPTER II 


REVIEW OF LITERATURE 

Crop surveys have been made by man since the raising of domestic 
crops was first started in the Middle East. Various methcdti of identi- 
fying and quantifying domestic crop production have been used in years 
past. Today, much of the world *s crop production is surveyed each 
growing season and many political and economic considerations are 
affected by the reports of probable crop production. Consequently, the 
more accurate the survey, the more accurate are the decisions which are 
based on crop production information. 

At present, a great amount of money is spent on crop surveys. 

Most surveys are of the "statistically random sample" type In which a 
small percentage of the cropland is checked by a field observer. Infer- 
ences are drawn from these observations about overall crop production. 

Multiapectral Reflectance from Crops and Soils 

In general, when electromagnetic radiation strikes a crop canopy 
or a soil surface, the following phenomena will occur: (1) a portion 

of the radiant energy may be directionally reflected (Kumar, 1972) but 
it ia most likely that it will be scattered reflectance; (2) a portion 
of the radiant energy may be absorbed and later emitted at a different 
energy level; and (3) a portion of the radiant energy may be .transmitted. 
In short, the total amount of radiant energy striking a crop canopy or 
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£t soil stirface is equal to the aniount reflected plus the amount absorbed 
plus the amouat trenaiaitted* Kumar (1972) has written an excellent 
review on reflectance from plants and soils end the following para- 
graphs are adapted from his discussion of reviewed literature. 

Interaction of Light with a Plant Canopy 

The analysis of reflectance from a plant canopy is extremely 
difficult because many vatiablea are involved* The moat important of 
these variables are: 

1. Absorption by oxygen, carbon dioxide and water vapor reduce 
incoming solar radiation in certain wavelength bands which 
decreases the accuracy of measurement* 

2. Illumination from the stan varies In intensity with numerous 
conditions* 

3* Kadiance from field plants Is affected by plant geometry, back- 
round soil reflectance, etc. 

4* The intensity of the sun has a maximum at about 0.5 vm, 
falling off rapidly at shorter and longer wavelengths. 

et. al. (1966) have shown that near infrared spectrophoto- 
meter studies of single leaves can be very misleading for predicting 
reflectance from crops. Near infrared light transmitted through the 
top of the crop canopy changes in light quality because of multiple 
internal reflections occurring within the canopy. Some radiation la 
scattered between leaves of a plant canopy by multiple reflection so 
that the reflectance (albedo) for the canopy as a whole is less than 
for single leaves and seldom exceeds Z5X» The amount of scattering 
increases with the irregularity of the leaf surface and with solar 
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elevation, because sunlight penetrates further into the canopy as the 
sun approaches the zenith. Davis (1957) has shown that the reflectance 
of grass varies with the altitude of the sun. His values varied from 
22% at noon to ab >ut 43% at sunrise and 48% at sunset. 

Existing theory of diffuse light propagation la not limited to one 
or two parameters. Silberstein (1937), for example, increased the 
nundier of parameters to three; an absorption coefficient and coefficients 
for both forward and backward scattering. In 1967, Duncan et. al, 
developed an elaborate theory for the penetration of direct ai.d diffuse 
sunlight through a foliage composed of many layers v f leaves with known 
orientation area, reflectance and transmittance characteristics. The 
controlled variables of the model ares leaf area, leaf angle, vertical 
position of layer of leaves, light reflected from leaves, light trans- 
mitted through leaves and the physiological relationship between 
illumination and photosynthesis. The variables of the environment are: 
elevation of sun above the horizon, solar intensity and sli^light 
brightness. Examples of computer simulations of hypothetical and real 
problems have been presented. Another elaborate model of a plant 
canopy was proposed by DeWlt (1965) . 

Anderson and Denmead (1969) have described the method for easy 
calculation of the flux densities of direct and diffuse radiation on 
inclined foliage in model plant stands. The stands are composed of 
randomly oriented, constantly Inclined, flat foliage surfaces. The 
calculations require knowledge of the flux densities of direct and 
diffuse radiation of a horizontal surface above the stand, foliage 
inclination, foliage area index and solar attitude. For direct 
radiation, the effects of changes in foliage inclination angle on the 
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average radiation received by the foliage are shown to depend strongly 
on solar altitude, and tine of the day. Relatively large differences 
can ejdlst between stands of different foliage inclination. There are 
only Btoall differences between surfaces of different inclination in 
the receipt of diffuse radiation^ particularly at the top of the stand. 

Allen et. al. (1970) have generalised and interpreted Duntley 
equations (1942) to account for the diurnal nature of near-infrared 
radiation measured in a com canopy. The Duntley optical coefficients 
associated with the specular component of light were assumed to vary as 
the secant of the aun*s zenith angle. Generalization of the Duntley 
relations was required in order to predict values of irradiance within 
the canopy and to account for the effect of background reflectance from 
the soil. Five independent measurements of canopy irradiance suffice 
to determine the Duntley parameters. Twenty-four measurements of trans- 
mittance within the canopy were used, however « to obtain a least squares 
calculation for the best fit of the Duntley equations to Irradiance 
within the com canopy. The Duntley equations fit the experimental 
results within a standard deviation of 3.2% for a period from noon to 
sundown. The best fit to near - Infrared transmittance measurements 
occurs when zero absorptance is assumed for the canopy. The Duntley 
equations reduce to a three-parameter representation for the special 
case of no absorptance. Other models of a leaf - Melamed Theory, and 
plate theory for a compact and a non-conpact leaf have not been applied 
to a plant canopy thus far. 



spectral Properties of Soil 

The pattern of spectral reflectance for soils is considerably 
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different from that for plants (Hyers and Allen, 1968) , In the near 
Infrared region, very substantial contrasts occur In reflectance 
betifeen different crop species and soil types. Spectral reflectance 
contrasts for soils, which show up as tone contrasts on photographs, 
are substantial throughout the spectral range of reflected solar energy 
(about 0.25 pm to 3 pm). 

Krlnov (1947) made the most extensive previous measurements of the 
reflectance of natural surfaces of soils, sands, and vegetation and 
showed that the reflectance of soils and sands increased monotonically 
with Increasing wavelength throughout the visible and near Infrared 
(to 0.9 pm). Bauer and Dutton (1962) observed the albedo values of 
agricultural areas, wooded hills, frozen lakes, and all of these areas 
covered by snow at semlregislar Intervals. The Instrumentation installa- 
tion consisted of two Eppley pyrhellometers and a Kipp and Zonen 
hemispherical solarimeter, mounted on a light aircraft, for measuring 
solar and hemispherical sky radiation, and reflected radiation. Values 
between 10% and 20% were observed over agricultural land areas in snow 
free seasons. With snow, the albedo values were as high as 80% over a 
frozen lake and as low as 50% over wooded hills. ‘ Gates (1964) has 
reported the spectral reflectance of some soil types. 

Several iJivestlgators have noted the so-called color effect on 
soil temperatures* The elevated daytime temperatures of dark-colored 
soils is attributed to their greater absorption and thus less reflectance 
of solar radiant energy. 



Several authors (Shockley e£ al. (1962), Kortum et al. (1963), 

Bowers end Hanks (1965), Orlov (1966)) have concluded that increasing 
particle diameter of soils results in a decrease of reflectivity. 

The conclusion is correct only for the laboratory case of dispersed 
soils. Zwennan and Andrews (1940), working with enameled surfaces, 
stated that spectral intensity of reflected radiation varies inversely 
with particle diameter. Orlov explained that the artificial breakdown 
of aggregates usually leads to an increase of the reflection coefficients 
caused by the character of the mutual position of aggregates. Fine 
particles fill the volume more completely and give a more even surface. 
Coarse aggregates, having an irregular shape, as a rule, form a very 
comples surface with a large number of interaggregate spaces (pores, 
cracks, etc.). Steiner and Gutermann (1966) described soil investiga- 
tions by Belonogova and Tolcbel*nikov (1959) and reported chat a 
decrease of grain size results in an Increase of reflectance, caused by 
greater scattering and lower extinction of light passing through the 
particles. Also, the area covered by microshadows occurring between 
particles under oblique illumination becomes smaller with decrease in 
grain size. They also demonstrated that the reflectance of soil miner- 
als depends on their dispersion in the soil. Structureless soils 
reflect 15% to 20% more light than soils with well defined structure. 
Reflectance varies with particle diameter but the shape of the spectral 
curve remains the same. Measurements by Coulson (1966) of the 
reflecting and polarizing properties of various soils, sands, and 
vegetation in the visible and near infrared spectral regions showed 
that dark surfaces polarise the reflected radiation strongly while 



highly reflecting surfaces have relatively weak polarizing properties* 

He found that reflectance of mineral surfaces increases with increasing 
angle of incidence and with Increasing wavelengths to about 2.2 pm* 
Atmospheric scattering affects principally the reflectance from short 
wavelengths and dark surfaces. Other factors such as locations acidity 
and past history of the soil cause difference in the multispectral 
responses received for any given species of vegetation grown on the 
soil* 

Shockley et al. (1962) reported the Influence of soil moisture and 
bulk density parameters on reflected energy In wavelength range 1.4 pm 
to 5.0 pm. They demonstrate the value of a soil moisture signature in 
identifying soils. Obukhov and Orlov (1964) stress that wetting and 
pulverization of the soil surface bring the reflectivity of soils 
closer to each other. Because of this, the most contrasting photo- 
graphs can be obtained at a low moisture content. A low contrast can 
also be expected with sun at a high angle above the horizon. Myers and 
Allen (1968) reported that wetting the soli in undisturbed and disturbed 
conditions substantially reduced the reflectance. 

Bowers and Hanks (1965) show that surface moisture content , organic 
matter, and particle size strongly Influence reflectance. Reflectance 
was found to decrease as moisture Increased. The staff of LARS at 
Purdue University obtained the spectral measurements of 250 soil 
samples. Ten different soil textures, four drainage profiles and three 
major soil horizons were represented in these samples. The mean 
spectral cprves for the clay soils at two different moisture levels and 
sandy soils at three different moisture levels are shown on page 84 of 
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LARS Vol, IV (1970) . These curves show a very large decrease in 
reflectance In the visible and near infrared range with Increase In 
iQolsture, Johannsen (1969) studied in greenhouse and field environ- 
ments the soil moisture - plant moisture relationships and the effect 
of these relationships on reflected and emitted energy from the soil 
and plant surfaces* 

f 

Crop Recognition Using Automatic Data Processing of 
Multiapectral Reflectance Imagery 
Most multlspectral reflectance imagery is obtained by an optical 
system measuring energy from discrete wavelength bands. The range of 
the energy measured is usually from about 0.4 micrometers (ym) to about 
2,6 mlcromsters (ym). This ta the visible and near infrared portion 
of the electromagnetic spectrum in which the reflected energy is most 
dominant. Both photographic and optical - scanning systems are used* 

In general^ the photographic systems are better adapted to photo- 
interpretation analysis and the optical - scanning systems producing 
digital output are better adapted to quantitative analysis by automatic 
data processing systems using computers. The system at LARS is of the 
latter type* 

In 1968, and again in 1970, the LAS^ staff reported that crops such 
as wheat, com, soybeans and hay could be accurately identified, using 
computer-aided processing of multiapectral reflectance imagery taken 
from an aircraft. Anuta and MacDonald (1971) reported on crop 
identification using digitized multiband satellite photography. Their 
results were only somewhat encouraging as they found that crops 
STdilblting a low amount of ground cover were indistinguishable from 


bare soils*. 
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Swain (1972) discusBea pattern recognition applied to digital 
imagery as a basis for remote sensing analysis. Hoffer (1973) discusses 
the use of automatic data processing (by computer) to analyze multi- 
spectral scanner data for land-use considerations* He concludes that 
automatic data processing is not only feasible but may become necessary 
when analysis Involves large spatial areas and temporal observations 
at the same time. Hall et al. (1974) have reported on problems 
encountered in attempting to use ERIS-1 imagery for crop identification 
studies. They found that location of field boundaries, field size and 
cloud cover caused the major problems* They concluded that, despite 
the problems encountered, analysis of EHTS-1 imagery by computer should 
be relatively cost effective* Landgrebe and the Z.ARS staff (1973) did 
an early evaluation of machine processing techniques of ERTS-1 data and 
discussed many different land-use situations in which ERTS-1 imagery 
might be useful* 

Bauer and Clpra (1973) reported that com and soybeans had been 
identified satisfactorily in northern Illinois using ERTS-1 multiband 
imagery* Computer processing of the multispectral imagery was used as 
the means of analysis. They obtained a overall accuracy for crop 
identification of 83%. They also found that using multitemporal imagery 
obtained during the growing season improved the identification of 
"other"; "other" being features other than corn or soybeans. The area 
estimates of com, soybeans and "other" obtained from the ERXS-1 
classification agreed closely with the USDA estimate. 

Williams et al. (1973) found that wheat could be identified by 
photointerpretation methods for a test site in Finney County, Kansas. 
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Using two classes, wheat and noa*^heat, and only the 0» 60-0. 70 ym wave- 
band of EBIS-1 satellite Imagery, they obtained an observed accuracy 


of wheat identification of 89%. 

» t 

Three separate investigations in California by Draeger (1973), 
Johnson and Coleman (1973), and Thomson (1973) have reported that crop 
identification is feasible if large fields are available for use in 
training and for testing accuracies. All of the three studies reported 
Identification accuracies of around 80% using photointerpretlve 
methods of ERTS-1 imagery. Horton and Heilman (1973) found that in 
South Dakota com and soybeans could be identified accurately, about 
90%, by machine processing all four wavebands of ERTS-1 imagery. 


CHASTER III 


METHODS AND OBSERVATIONS 


Site DescriptionB 
Greeley County, Kan see 

Greeley County Is one of the western Kansas counties on the 
Colorado border and is centrally located between Nebraska and Oklahoma* 
Its area is about ?.00,000 hetares* It occupies part of the nearly 
level to gently rolling high plains region between the Arkansas River 
to the south and the Satokey Hill River on the north. 

Agriculture is the only Industry of Greeley County. Wheat and 
cattle are the main sources of income* Most of the county is cultivated 
and most of the farming is dry-land. There is a email amount of irriga- 
tion from deep-well sources. Most of the areas remaining in native 
grassland are found on slopes adjacent to natural drainage ways. 

The soils of the county are quite uniform in color and surface 
texture. They are quite dark when moist and light-gray when dry. 

About 93Z of the soils have a silt loam surface texture. The topography 
is gently undulating and the soils are quite susceptible to wind erosion. 

One striking cultural feature of the county is field patterns* 

Much of the land in wheat production is sectioned into long, narrow 
fields with the long aides of fields lying in an east-west direction. 

The prevailing winds in western Kansas are of a north-south direction. 
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By locating the fields perpendicular to the dominant wind direction, 
the wheat and wheat stubble acts as a natural barrier' or windbreak. In 
this manner, much of the wind erosion can be retarded. 

Hard red winter wheat is grown in Greeley County. It la planted 
in the fall, germinates, and tillers before winter, goes through 
hardening and vernalization in the winter, and has its major growth, 
flowering and maturity in the spring. Wheat harvest in the county 
usually starts about the third week of June, Hard red winter wheat is 
grown in the county because the annual rainfall is only about 40 cm per 
year. Fallowing on about one-half of the wheat land is practiced. When 
land is "fallow” it is allowed to lay idle for one year while it 
accumulates moisture in the soil and has the weeds controlled. The 
following year, the “fallow" land is planted in wheat and the land 
which was in wheat the previous year is "fallow". 

Lubbock County, Texas 

Lubbock County lies in the Southern High Plains in west-central 
Texas. It is slightly south of the Texas Panhandle and is near the New 
Mexico border. The county encompasses an area of approximately 230,000 
hectares. Yellow House Draw bisects the county on a diagonal from the 
north-west to the south-east comers. Most of the topography in areas 
adjacent to the draw is quite rolling. The remainder of the county la 
nearly level to gently rolling xrtth the most pronounced topographic 
features being playas} playaa are small depressional areas in the 
landscape which catch and hold the runoff water for a short period of 
time. After a rain the play as resemble small lakes. A laka, in 
Spanish, Is playa or beach. 



The main agricultural products of the county are cotton^ grain 
sorghum and cattle. Almost all of the cotton and grain sorghum la 
irrigated with water from deep wells. Small amounts of dryland wheat 
are grown. A few soybeans are *'“so grown. Much of the land near the 
draws and dralnageways is in native grasses which are used to pasture 
cattle. 

The climate in the Lubhoclc area it? also semi-arid. The water 
loss due to evaporation and transpiration make it necessary to Irrigate 
most summer crops. The soils are quite heterogeneous and esdiibit some 
limitations to production. Salt accumulation, low natural fertility 
and periodic wetness (In playas) are some of the soil deficiencies. 

The main agricultural crops in Lubbock County are "short-day** 
plants, flowering after the days in summer start to get shorter. Grain 
sorghum is planted about the middle of June. Cotton is planted near 
the end of June following grain sorghum planting. Almost all of these 
crops are planted on ridges. Irrigation water can then flow down the 
furrows between the ridges to Irrigate the crops. Both crops are 
harvested In late fall. 

One observation that can be made in an area which uses furrow 
irrigntlon practices is that fields may be irregular in shape and size. 
The land must be leveled and provided with a slight, unidirectional, 
cvonstant slope to Insure that the irrigation water will reach all of a 
field uniformly. Fields are leveled in a manner which allows the least 
amount of soil to be moved. Therefore, original topography is taken 
into consideration and the resulting fields may have irregular shapes 


and different sizes. 
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Acqulaltion of Multlspectral Data 
The sniltispectral data for this study was collected by ERTS-l 
which was launched Into a polar orbit In July of 1972. The orbit is at 
an angle of .14*^ with the longitude of the earth. This 14* angle 
provides a sun-synchronous orbit which places the satellite over Its 
target at approxitaately 9:30 A.K. local time. 

The multispectral scanner on ERTS-1 measures reflected energy in 
four wavelength bands. The bands are as foUowe (ERTS-1 Data Users 
Handbook, 1972): 

0.50 - 0.60 micrometers (pm) green 

0.60 » 0.70 micrometers (pm) red 

0,70 - 0.80 micrometer 8 (pm) infrared 

0.80 - 1.10 nicrometera (pm) infrared 

The multispectral scanner on ERTS -1 is an optical-mechanical scanner 
with a field of view of approximately 185 kilometers. The images or 
scenes obtained cover about 185 by 185 kilometers or a little over 
3,400,000 hectares. 

The multispectral data for this study was obtained on computer 
compatable digital tapes from the Goddard Space Flight Center at 
Greenbelt, Maryland, The data was received at LARS as part of the 
multispectral data for NASA Contract NAS5-21785. Only one scene was 
analysed for the test site in Kansas. Four scenes covering the growing 
periods of cotton and grain sorghum were obtained for the Texas test 
site. 

Table I lists the ERTS-1 Images used. All of the multispectral 
data was obtained during the summer and fall of 1973 by the ERTS-1 
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Table 1« ERTS-1 Images used for crop identification. 

Location Pate ERTS"1 Serene IP 

Greeley County, Kansas June 19, 1973 1331 - 16571 

Lubbock County, Texas June 18, 1973 1330 — 16531 

Lubbock County, Texas July 24, 1973 1366 - 16521 

Lubbock County, Texas August 11, 1973 1384 — 16524 

Lubbock County, Texas October 22, 1973 1456 - 16523 

satellite. Date was collected on June 19, 1973 for Greeley County, 
Kansas and on June 18, July 24, August 11 and October 22, 1973 for 
Lubbock County, Texas. Each image contains about 7.7 million data 
points. This means that the 'response represented by one data point in 
each waveband is the integrated response from 0.44 hectares. Only Che 
portion of each image covering the respective test areas was analyzed. 

Ground Control Information 

Ground control Informatio- is that data used to train a photo- 
interpreter or a computer how to recognize certain characteristics 
(spatial or spectral) of a scene as having a specific feature. For 
example, if an analyst wishes to Identify com fields in a heterogeneous 
agricultural scene photographed from an airplane, a knowledge that some 
specific fields located in the scene contained com can greatly aid In 
identifying correctly the remainder of the com fiel.ds. For these 
purposes, some ground control information was essential for this study. 

Ground control data for Greeley County, Kansas was from three 
sources. One was simply the Statistical Heporting Service (USDA/SRS) 
estimate of the total amo\mt of wheat grown In Greeley County in 1973. 

A second source was in situ observations of some selected fields by the 
Cooperative Extension agent in Greeley County. About 25 fields around 
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Tribune* Ranaae were observed by the Ckioperative Extension Agent. He - 
reported whether fields had been In wheat* fallow* or In permanent 
pasture for the 1973 crop year. Both this data and the SRS estimate 
of wheat were used to verify the accuracy of the Identification of 
wheat In Greeley County. 

The third source of ground control data and the only ground 
control used for training was aerial color infrared photography. This 
photography was obtained by an aircraft flying at an altitude of 9500 
meters on a south to north flightline over Kansas State Highway 27* 
which passes through Tribune* Kansas and is near the middle of the 
county. The photography was taken on May 14* 1973 when the wheat was 
a lush green* the permanent pasture was starting to become green and 
the fallow land was bare. 

Interpretation of a scene photographed with color infrared film 
requires a knowledge of some basic characteristics of the film. The 
film is sensitive to a near infrared waveband (0.72 - 0,92 micrometers) 
and requires a filter on the optics which screens out the blue wave- 
length radiation. Normal color film depicts a blue river with a blue 
color and a green tree as a green color. Color Infrared film displays 
blue targets with a black color* a green* non-living target as a blue 
color* a red target as a green color* and a lush green vegetative cover 
as a red color. The red color on color Infrared film corresponds to 
the 0.72 - 0.92 micrometer waveband. Compared with other types of 
cover* green plants reflect more highly the near Infrared energy, thus 
the usefulness of color Infrared photography when working with green 
vegetation Is apparent. 
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Ground control data for Lubbock County, Texas was obtained in situ 
by farmers In the area and by M. F. Baumgardner and J. An Henderson 
for LARS. Most of the data was on crop type and crop conditions. Soce 
low altitude aerial photography was taken by Henderson. 

Four farmers collected data which was used in this study. Figure 
1 indicates the locations of farmer ground control information sites 
with arrows. Each farmer observed all of the fields on both sides of a 
paved road near his home. The length of the segment he surveyed for 
each ERTS-l pass was ten to twelve kilometers. Each farmer reported 
land use (crop), planting pattern, growing conditions, percent ground 
cover, stage of crop residue, soil conditions, crop conditions and row 
direction for each field in his segment. The voluntary and diligent 
effort of these four farmers was greatly appreciated. Other ground 
observers did not collect ccmiplete information and their observations 
could not be used. 

During the first week of July In 1973, Baumgardner and Henderson 
took low altitude aerial photography of thirty-six (36) road Inter- 
sections in Lubbock County (See Figure 1) » These road intersections 
were located on three lines running north to south through the county; 
each line was along a north-south county road and each line covered 
about one- third of the county. The urban area near the city of Lubbock 
was avoided. Each intersection was marked on a county road map so that 
the Intersection^ could be located when on the ground. 

The day after the aerial photography was collected each of the 
intersections was visited. The land use and other features like those 
noted by the four farmers were recorded for each . the four corner 
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Figure 1, Location of ground information used in the analysis of Lubbock County, Tesas. 
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fields at each InterBectlon. This gave an additional 144 fields of 
ground control to be used in training or to test classification 
accuracies* This gave a total of about 300 fields with ground control 
information in the Lubbock teat area^ though many of the fields were 
too email for effective use* 

For a more complete description of the method of obtaining ground 
control for this experiment the reader should refer to the final report 
for tIASA Contract NAS5-21785 by Baumgairdaer (1974) » 

Analysis of Multiapectral Satellite Imagery 

The LAKSYS software system was used to analyze the multlspectral 
satellite imagery. LARSY3 is a package of computer programs which has 
been designed to analyze and display remotely sensed multlspectral data. 
The use of these programs is discussed by Phillips (1974) and Hoffer 
(1973), The computer used was an IBM 360-67. 

Eight separate processing algorithms were used in this study: 

(1) GEMCOR, (2) IMAGE0ISPLAY, (3) CLUSTER, (4) STATISTICS, (5) 
SEPARABILITY, (6) CLASSIFYPOINTS p (7) PRINTRESULTS, and (8) PHOTO. The 
first five algorithms were used to analyze data from both test sites. 

The last three were used only for the Greeley County test site. The 
results of the Lubbock County analysis were such that analysis beyond 
SEPARABILITY was not necessary. 

The first step in the analysis process was to correct the multi- 
spectral data for geometric distortion. Due to the multlspectral 
scanner geometry and the heading of ERTS-1, considerable spatial dis- 
tortion occurred. The algorithm GEMCOR, reported by Anuta (1973), was 
used to correct the distortion and adjust the scale of the multlspectral 
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data* By locating tha approximate comer coordinates of the counties 
in all of the scenes and supplying them to GEMCORp the area ^^thin the 
designated coordinatea was corrected for geometric distortion and scale. 
Approximate cdmer coordinates ffere located using 70 millimeter trans- 
parencies of the 0.60 - 0.70 pm waveband of the multlspectral imagery, 
provided with the digital Imagery. The corrected sets of digital 
imagery were used fo^ all further analysis. 

Next, the corrected digital imagery was displayed, one scene at a 
time. The Digital Image Display System displays the image on a black- 
and-white television screen. An interactive capability to edit, 
annotate or modify the image is provided through a light pen and a 
program function keyboard. An additional photographic copying capa- 
bility is also available* 

The computer program which allows the interaction of the Digital 
Image Display System and the data set Is IMAGEDISPZJVY. The data in 
each Individual waveband is partitioned into 16 levels and these levels 
are displayed on the screen as gray levels, low values being dark and 
high values being bright. The program also provides many other functions 
such as outlining fields with the light pen and obtaining the coordin- 
ates in the multlspectral data, magnifying the image on the screen, and 
many other features. 

The corrected scenes were displayed and the ground control sites 
were located In the multlspectral Image, usually using the 0.60 - 0.70 
pm and 0.80 - 1,10 pm wavebands. A rather large area around a ground 
control site was outlined since exact location of a ground control site 
on the display system is most difficult at times, especially when 
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dealing with agricultural crops. For example, in Greeley County an 
area of about 8 km by 32 km was outlined and coordinates were obtained. 
Also, in Greeley County, the exact corner coordinates of the county 
were located in the oultlapectral Imagery and t^JOse coordinates 
recorded. In Lubbock County, a aimilar procedure was used in outlining 
the areas around the ground control sites. 

After the portion of tsultispectral data taken over the ground 
control sites was located, the CLUSTER algorithm was used to produce a 
Ttap (computer printout) of the area. CLUSTER is an unsupervised 
classifier, a pattern recognition tool that groups data vectors into 
an arbitrary number of spectrally distinct classes. To enhance the 
field boundaries which were somewhat indistinct and undlstinguishable 
on the display system, the data was clustered into eight spectral 
classes. Each data point within the ground control site was assigned 
to one of the eight spectral classes by CLUSTER and was displayed on 
the computer printout as one of eight different S 3 ^bols. 

Unique features such as odd-shaped fields, lakes and road inter- 
sections were used to match the CLUSTER map with the ground control 
data. Then definite fields with known crop type could be located 
spatially in the CLUSTER map and field coordinates were obtained for 
each field selected. County road maps, ground control field maps and 
low altitude airphotos were all useful for precise location of fields 
in Lubbock County* In Greeley County, the field size and shape could 
be seen in the color infrared aerial photography. 

The land use or crop type in Greeley County was also determined 
from Che color infx >red photography. On Hay 14, 1973 the wheat fields 


In the color infrared photography were bright red, pasture fields were 
light pinh, and fallow land was black or greenish-brown. 

Temporal overlay capability, the aligning of data sets of the same 
area differing only in time, was employed for the Liibbock County teat 
site. This alignment, or spatial registration of multitemporal data, 
matched the coordinates of a given point on the ground for all four 
scenes of data used. This made locating ground control fields in 
Lxibbock County necessary only once and eliminated some possible sampling 
error. 

The STATISTICS processor was used to obtain mean vectors and covar- 
iance matrices for the different classes of crops or land-use selected 
from each test site. A class mean vector was calculated by averaging 
the response from all of the data points within all of the fields used 
for training for a specific class (crop or land-use) . All four wave- 
bands were considered. The relationship between the waveband responses 
for a specific class was shown in the covariance matrix. 

Five classes of land-use were selected for training in Greeley 
County. They were wheat, pasture and three types of fallow land. 

Three classes of fallow land were selected because of differences in 
cultivation. Some fields were recently cultivated, others were not 
cultivated when the mxxltispectral imagery was obtained, and some were 
weedy. 

In Lubbock Coimty, Texas a number of classes were defined but 
cotton and sorghum were the two classes of interest. The main objective 
\t this test site was to identify cotton and sorghum and differentiate 
the two crops. For comparison purposes classes of permanent pasture. 
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temporary pasture^ water^ and urban cosmtuniCies were chosen* A 
separate set of mean vectors and covariance matrices was calculated for 
each of the EETS*>1 images. 

The vectors and covariance matrices for each scene were used 
as input for the SEPARABILITY processor, an algorithm which measures 
the statistical distance or separability of the class vectors. The 
processor considers a specified set of wavebands of data (in this study, 
all four wavebands of the ERTS-1 imagery were specified) and computes 
a transformed divergence value (Swain, 1973) for all possible combina- 
tions of classes. It has been essperimentally observed that a minimum 
value of transformed divergence of 1600 is required if classes are to 
be considered separable. Values lower than 1600 tend to indicate that 
the two classes being considered are similar and the probability of 
discriminating between them accurately is quite low. Values higher 
than 1600 Indicate the classes are separable. The maximum value 
possible is 2000. 

The SEPARABILITY processor is a good method to check to see If the 
training classes that have been selected will produce acceptable 
classification results. Classes which are deemed Inaeparable by 
SEPARABILITY are not likely to produce accurate classification when 
used in the CLASSIFYPOIKTS algorithm. 

The CLASSIFYPOINTS processor uses the mean vectors and covariance 
matrices for training classes calculated by STATISTICS to perform a 
maximum likelihood Gaussian classification, data point by data point, 
for a specified set of data. The classifier (pattern recognition 
algorithm) compares the response at each data point specified with the 
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statistics of the traiaing classes and makes a maximum likelihood 
estimate as to which class it belongs. In Greeley County » the data set 
specified was the area in the data which was within the county boundrias. 
No classification was perfom^d for Lubbock County because the results 
of SE^ABABILITY Indicated that the resulting classification would be 
highly inaccurate. 

The ?RXNTHES1ILTS processor has two main functions s (1) to display 
the polnt-by-point classification of the specified data as an 
alphanumeric map-like printout on a lin«i printer, in which the user 
selects the symbols to be used for each of the different classes, such 
as C for cotton or for wheat | and (2) to produce a quantitative 
evaluation of a classification in tabular form. A table listing how 
many data points fell into each specific class is one product. Also, 
coordinates of fields of knom land-use or crop type (not fields used 
in training) can be input into the program as test fields. PRINTRESULTS 
will compare the classified points within the test fields with their 
known class and will compute a table of classification accuracies. 

This is the most T?idely used method, for crop studies, to evaluate the 
accuracy of the CLASS IFYPOINTS procedure. 

Another method of displaying the results of the CLASSIFYPOIHTS 
procedure is the photographic capability of the Digital Image Display 
System, A program called PHOTO causes the different classes of the 
classification to be displayed at selected intensities on the Image 
display screen* A black and white photograph may be taken at this tlnie 
from the photographic image display screen. Color photography is also 
possible with PHOTO. The color for each class is selected from a color 
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chart and is specif led in the program. Then, by using a combination 
of filters and ordinary color film, color-coded classifications can be 
produced. 
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CHAPTER IV 

RESULTS AND DISCUSSION 


Greeley County. Kansas 

The mean and standard deviation of the relative spectral responses 
for wheat and other land-use types in Greeley County on June 19, 1973 
are listed in Table 2. All relative response values arc rounded to the 
nearest whole number. 

Table 2, Response of wheat and other land-use types in Greeley Cotmty 
(X and s). 


Land-Use Wavelength Band (pm) 


Type 










«*■ 

X 

8 

X 

s 

X 

s 

X 

8 

Wheat 

37 

3 

40 

5 

51 

3 

28 

2 

Pasture 

39 

2 

38 

2 

53 

3 

28 

2 

Fallow 1 

43 

2 

50 

2 

50 

2 

24 

1 

Fallow 2 

61 

4 

73 

6 

73 

4 

35 

2 

Fallow 3 

53 

3 

63 

4 

63 

3 

31 

1 


The means are plotted in Figure 2. The relative reflectance of 
fallow land is higher than wheat or pasture in the visible wavebands. 
The response in the 0,50-0.60 pm waveband is higher for pasture than 
for wheat. The same is true for the 0.70-0.80 pm Infrared waveband; 
however, the reverse is true in the other two wavebands: the response 

for wheat is higher than pasture. It would not be correct to compare 
the taeans from waveband to waveband because the wavebands are all 




Figtire 2. HuXtispectral response of wheat end other land-use classes in 
Greeley County. 
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calibrated individually. The response inversions for wheat and pasture 
and their large differences in response with the fallow classes 
suggested that the five classes were separable classes. 

Ihe standard deviations listed in Table 2 were considered rela- 
tively aiBall. This suggests that the classes selected have a small 
response variance in each waveband and the probability density func- 
tions for the classes would be less likely to overlap. Fallow 1 was 
believed to be freshly cultivated land. Fallow 2 was considered uncul- 
tivated, and Fallow 3 was thought to be uncultivated, weedy fields. 

The SEPAHABILITY processor affirmed the conclusion that the wheat 
could probably be classified (Identified) correctly. Table 3 lists the 
results obtained from SEPARABILITY indicating that the multispectral 
response of wheat was different enough so that there should be little 
confusion with the other classes considered. The average transformed 
divergence for all class pairs was 1921 (2000 is the maximum obtainable 
value) . The minimum divergence between wheat and any of the other 
classes was 1739, On this basis, the wheat should be identified 
correctly as wheat by the CLAS SIFYPOINTS processor. 


Table 3. SEPARABILITY results listing divergence of the five classes 
identified in Greeley County, Kansas using all four ERTS-1 
wavebands. 


Class Combination 
Compared 

Transformed 

Divergence 

Wheat vs. Pasture 

1739 

Wheat vs. Fallow 1 

1999 

1['7heat vs. Fallow 2 

1996 

Wheat V8, Fallow 3 

1995 

Pasture va. Fallow 1 

2000 

Pasture vs. Fallow 2 

2000 

Pasture vs. Fallow 3 

2000 

Fallow 1 vs. Fallow 2 

2000 

Fallow 1 vs. Fallow 3 

1995 

Fallow 2 vs. Fallow 3 

1487 
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When the entire county was classified, the classification was 
tested for correct identification of wheat and the other classes by 
introducing test fields of known land-use type. The test fields were 
chosen from the under flight color Infrared photography and from a 
limited number of fields visited by the local county agent. All test 
fields for fallow land were combined for purposes of testing the 
classification accuracy (i.e., no attempt was made ...). The percent 
correct classification was computed by taking the total number of data 
points classified correctly within the test field coordinates for each 
class and dividing that number by the total number of data points for 
that class. Table 4 lists the classification accuracies. 

Table 4. Classification accuracy for wheat, pasture and fallow In 
Greeley County. 


Class 

Name 

No. of samples 
per class 

% Correct 
Classification 

No. of samples classified with 
Wheat Pasture Fallow 

Wheat 

400 

97.0 

388 

4 

8 

Pasture 

318 

96,1 

9 

305 

4 

Fallow 

431 

97.9 

7 

2 

422 


These classification accuracies for tested fields were considered 
excellent. A further test was made of the classification. Multiplying 
the number of data points in the county classified as wheat by 0.44 
yields the number of hectares of wheat identified in Greeley County. 
Table 5 compares the results of a Statistical Reporting Service (SRS) 
estimate of wheat In Greeley County for 1973 with the number of hectares 
of wheat Identified by the analysis of the Imagery. The SRS estimate 
was obtained from the county extension agent. 
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TabXe 5* 1973 area estimates of wheat, fallow and row- 

crop, end pasture In Greeley County. 

Fallow end 

Wheat Pasture Row-Crop 
Source (Ha) (Ha) (Ha) 

SRS (USDA) 77,000 35,000 88,000 

Identification via ERTS-1 78,000 32,000 92,000 

The USDA/SRS estimate for wheat is only about 1% smaller than the 
amount of wheat Identified using ERTS-1 imagery. This close agr effluent 
of statistics carries through in the pasture class, the difference 
being about 4%. SRS row-crop estimates were combined with estimates 
of fallow (12,000 Ha and 76,000 Ha, respectively) and were displayed 
as one class in Table 5, At the time the ERTS-1 Imagery was obtained 
for Greeley County (June 19, 1973) any row-crops such as com or 
sorghum would have low ground cover and would look like the bare fallow 
fields. Identification of wheat and other land-uses by satellite 
imagery is similar to the estimate obtained using present SRS techni- 
ques for Greeley County, Kansas. It is impossible to determine with 
available data which of the two estimates la the better estimate of the 
actual situation. 

A visual aid or map is useful for observations of spatial distri- 
bution and cultural patterns of wheat in Greeley County. By using the 
PHOTO processor, a photographic map of the wheat classified in Greeley 
County was produced and is shown in Figure 3. The long, narrow fields 
of alternating fallow and wheat are easily observed. This cultural 
practice is used to control wind erosion. The dominant pasture areas 
appear along the drainage vays as they should. The area displayed in 
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Figure 3. Photographic map of wheat classification for 

Greeley County; color code is: white = wheat, 

gray “ pasture, and black » fallow land. 
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the piiotograph is about 200,000 hectares. The classification appears 
spatially reasonable as well and Illustrates the homogeneity of the 
agricultural practices in Greeley County. Under similar conditions, 
wheat, and perhaps other crop types, can be ^txpected to be identified 
correctly using satellite imagery. 

Regression analysis was used to help exp?^ln the reflectance 
character is tics of the various cover types snd the relationships among 
the several variables. A multiple regression model using the four 
wavelength bands as independent variables was adapted from Draper and 
Smith (1966). The dependent variable, land-use or surface condition, 
was coded with a five factor, orthogonal polynomial as follows: 
fallow 3“ -2, fallow 2« -1, fallow 1“ 0, wheat ■ 1, and pasture ■ 2. 
The practice of coding qualitative variables with orthogonal poly- 
nomials reduces bias In the regression model (Anderson, 1974) , The 
regression model used was: 

Y - b^ + bi + b^ Xg + bs X3 + b^ + e where: 

Y is the dependent variable for land use; 

is the intercept when X^^ * X2 " - 0; 

^1* ^2* ^3* ^4 coefficients in the regression model; 

is the response in the 0,50 - 0.60 pm waveband; 

X^ is the response in the 0,60 - 0.70 pm waveband; 

is the response in the 0.70 - 0,80 pm waveband; 

X^ is the response in the 0.80 - 1.10 pm waveband; 

and e is the random error tern distributed normally with mean 

2 

equal to zero and variance equal to o . 
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A plot of land-use with the response from the 0.60-0.70 ym wave- 
length band Is shown In Figure 4. This simple linear regression for 

2 

one waveband yielded an r ™ 0.757 which is very good for one waveband* 
This says that about 767S of variability In spectral response from the 
land-use classes Is explained by the response noted with the 0.60 - 
0.70 ym wavelength band. 

A prediction equation using multiple regression was formed In a 
stepwise manner* A step-^lee regression procedure used enters the 
Independent variable which explains the most variability into the 
equation first, then adds the second beat Independent variable, and 
continues to add Independent variables to the prediction equation until 
It has no more to add or until the next independent variable makes no 
significant contribution to the model (when the error sum of squares 
Is not reduced significantly) . The prediction equation formed for 
land-use conditions in Greeley County was formed without the 0.70-0,80 
ym waveband as it contributed nothing to the model: 

Y 1.26 + O.lSXj^ - 0 . 23 X 2 

2 2 
The multiple r for this regression model was r *0.85, Thus the 

three variable model explains 10% more than the linear model in Figure 

4 and should be quite useful In predicting what land-use type occurs 

if the spectral response is known* By substituting the spectral 

response for the respective independent variables, the type of surface 

condition can be computed and can be expected to be accurate about 83% 

of the time. 



CODED 

- 2-10 1 2 

SURFACE CONDITION 


Figure 4* Plot of surface conditions and the spectral response 
from the 0,60 - 0,70 m wavelength band. 
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Lubbock County, Texaa 

In this test site, cotton and sorghum were the tvo crops that were 
Investigated. The nultlspectral response measured by ERTS-1 was 
obtained from a set of training fields for each of the four ERTS-1 
images obtained during the growing season. Table 6 lists the results 
for cotton and sorghum by image date and by wavelength band. The means 
and standard deviations (s) are products of STATISTICS and are 
rounded to the nearest integer value. 

Table 6. Means (X) and standard deviations (a) for cotton 
and sorghum by image date and waveband. 


Date 

Waveband 

vm 


Cotton 
X s 

Sorghum 
X 8 

j 

0.50-0.60 

1 

38 

4 

37 

3 

U 

0,60-0.70 

2 

48 

6 

46 

6 

N 

E 

0.70-0.80 

3 

53 

5 

51 

5 


0.80-1.10 

4 

27 

3 

25 

3 

J 

0,50-0.60 

1 

42 

5 

41 

10 

U 

0,60—0.70 

2 

48 

10 

43 

13 

L 

Y 

0,70-0.80 

3 

68 

5 

64 

9 


0.80-1.10 

4 

34 

3 

32 

4 

A 

U 

0,50-0.60 

1 

32 

3 

32 

4 

6 

0,60—0,70 

2 

32 

8 

33 

9 

U 

s 

0.70-0.80 

3 

59 

6 

59 

6 

T 

0.80-1.10 

4 

32 

4 

32 

4 

0 

c 

0.50-0.60 

1 

28 

3 

28 

4 

T 

0.60-0.70 

2 

33 

6 

33 

7 

0 

B 

0.70-0.80 

3 

41 

5 

39 

6 

E 

0.80-1.10 

4 

21 

2 

21 

3 


R 



A suspicion developed immediately upon observing Che means and 
standard deviations of Table 6 that perhaps cotton and aorghum could 
not be differentiated under the conditions that existed in Lubbock 
County. The modified bar graph ahown as Figure 5 illustrates the close- 
ness of the means and the overlap of the data distributions. The 
means are quite close in most cases and are the same in soma wavebands. 
Furthermore, the overlap suggests that It is improbable that cotton 
can be separated spectrally from sorghum in this study. 

The SEPABABILITY processor was used to quantify the difference in 
the multlspectral response of cotton and aorghum. Values of 300 to 
400 were obtained as the measure of separability of cotton and sorghum 
for Individual dates and for combinations of dates by SEPARABILITY. 

An acceptable value for class separation is on the order of 1600 or 
greater so the cotton and sorghum were considered quite inseparable, 
spectrally. Also, the cotton and sorghum were found to be inseparable 
spectrally from all other cultivated crops in the Lubbock test site. 
Therefore, no classification was made of Lubbock County as it was 
unlikely that accurate identification of cotton and sorghum could be 
achieved. 

A more conventional means of determining the closeness of the mean 
spectral response for cotton and sorghum is Analysis of Variance. The 
model for the analysis is a split-plot design (Anderson, 1974). The 
model is: 

’'ijk - M + + Cj + 4. + CB^^ + DCB^^ 

+ where: 



RELATIVE RESPONSE 



WVEBANDS (l=0.5-0.6/im, 2=0.6-0.7/Am, 3=0.7-08/im, 4=0.S-I.I^) 


Figure 5. Mean response and standard deviation for cotton and .sorghum from their 
respective ERTS-1 images by waveband. 
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^Ijk 


M 

D. 


<i) 


''J 

DC 


ij 


D& 


ik 


CB 


Jk 


Is the response of the kth wavelength band from the Jth 
crop type within the ith scene of data; 

Is the overall mean 5 

Is the variation due to the ith scene of data; 
is restriction error due to each scene being a single 
set and not replicated; 

is the variation due to the jth crop type; 
is the variation due to the Interaction of the ith scene 
with the jth crop; 

is the variation due to the kth wavelength bond; 
is the variation due to the interaction of the ith scene 
with the kth wavelength band; 

is the variation due to the Interaction of the jth crop 
type with the kth wavelength band; 

is the variation due to the interaction of the ith scene 
with the jth crop type and the kth wavelength band; and 
is the error. 

As there Is only one observa- 
tion per cell, the and the terms cannot be separated and 

will be called error in the results of the analysis of variance in 
Table 7. The F-test for significant difference gives no significant 
difference for the variation due to crop type C^j)> This is computed 
by dividing the mean square of by the mean square for This 

gives a calculated F-value (12.5/4.1) of 3.05 with one and three (1,3) 
degrees of freedom. The tabular F-value for one and three degrees of 
freedom and at alpha ■■ 0.05 (F^ ^ 


DGB 


ijk 


^(ijk) 

The means shown in Table 6 are the Y 


is 10.13. Since the calculated 
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F-vaXue is eisaller than the tabular F-value, the variation in spectral 
response due to crop type is determined to be not significant at a 952 
level of confidence. The mean square for had to be used as the 
denominator in the F-test Instead of the mean square for error because 
the EKTS-l scenes must be considered a random variable. 

Table 7. Results of analysis of variance for cotton and sorghum. 


Source of 
Variation 

Degrees of 
Freedom 

Sum of 
Squares 

Hean 

Square 

Significance 
of F-test 


3 

1048.8 

349.6 

no test 


0 

0.0 

0.0 



1 

12.5 

12.5 

no 

“'ij 

3 

12.3 

4.1 

no 


3 

2974.5 

991.5 

yes 

“ik 

9 

456.3 

50.7 

yes 


3 

2.5 

0.8 

no 

error 

9 

4.8 

0.5 


TOTAL 

31 

4511.5 




All of the results from the Lubbock County test site indicate 
that cotton and sorghum cannot be identified or delineated using ERTS~1 
multlspectral imagery and pattern recognition techniques. This Is in 
definite contrast with the results obtained for the Greeley County 
test site. The explanation may be that the Lubbock area is hetero- 
geneous agriculturally. The fields are not uniform in efdior size or 
shape. There is also a broad variety of crop types. Most of the 
soils have a predominantly reddish-brown color which may dominate the 
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response characteristics of the crops when the ground cover Is low. 
Also* cotton and sorghum are usually planted about the sene time 
and have the same planting practices. Finally* although cotton and 
sorghum differ morphologically* they have about the same ground 
coverage while growing and the results suggest that the morphological 
differences are not great enough to be distinguished from satellite 
altitudes. The conclusion must be that cotton and sorghum cannot be 
identified accurately under conditions like those which existed at the 
Lubbock County test site using current ERTS-1 imagery. 
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CHAPTER V 
SUMMARY 

The oveirall objjective of the study was to identify important crop 
types In a semi-arid climate, using ERTS-1 muXtispectral imagery and 
pattern recognition techniques* Greeley County In western Kansas and 
Lubbock County in the western high-plains region of Teicas were the teat 
sites chosen* Wheat, cotton and sorghum were the major crops. 

Excellent results were obtained in Greeley County with wheat being 

identified correctly 97% of the time. The estimate of wheat for the 

whole county was within 5% of the USDA/SRS estimate* A multivariate 

regression equation formed for cover types in the county yielded an 
2 

r ■0*85. For a situation similar to that in Greeley County, ERT3-1 
multiband imagery appears to be a good tool for identification and 
area estimation of wheat provided that the multispectral Imagery la 
obzalned in late spring* 

The results for the more heterogeneous area of Lubbock County did 
not confirm the hypothesis that crops can be identified using ERTS-1 
multispectral imagery. In this study, cotton and sorghum could rot be 
separated from other cultivated crops using ERTS-1 Imagery collected 
during a complete growing season using Image dates taken in June, July, 
August and October* Heither were they spectrally dissimilar enough to 
be identified as two separate crops. These results are important 
because they Indicate that present EBTS-1 imagery may not be suitable 
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for ideatifylog crops ia areas with characteristics similar to Lubbock 
County. 

The results demonstrate the feasibility of using this space-age 
technology for obtaining crop production information if the crops are 
spectrally separable. Areas having uniform soils and cropping patterns « 
and relatively fe^ crops are most likely to meet these conditions* 

Soma knowledge of the land and its cropping patterns is essential for 
any crop survey, especially if current ERTS-1 laujgery is being used. 
Hultispectral satellite imagery currently available does not appear 
acceptable for identifying cotton and sorghum in areas with extremely 
mixed cropping patterns and soils. If, however, a sensor with greater 
spectral and spatial resolution, wavelength bands in the middle and 
thermal Infrared, end greater signal to noise ratio were available, 
it might be posaible to accurately classify crops under these more 
difficult situations. 


a 
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